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Note to the reader 



This paper is an alternative version of the pubhshed paper [Bu3], containing more 
results obtained in [Bui] and no numerical results. The paper is dated by the last 
modification of the file. 

1 Introduction 

This paper is devoted to the study of quasi-static brittle crack evolution. We 
work under the following assumptions: a linear elastic body, with or without initial 
cracks inside, evolves in a quasi-static manner under an imposed path of boundary 
displacements. During its evolution cracks with unprescribed geometry may appear 
and/or grow. 

The difficulty of brittle crack propagation problems consists in the nature of the 
main unknown: the crack itself, at various moments in time. The research in this 
field concerns mainly the constitutive behaviour of a brittle material, like the basic 
paper of Griffith [G]. The essential was stated in papers like: Eshelby [Es], Irwin [I], 
Gurtin [Gul], [Gu2]. 

In almost all the studies the geometry of the crack is prescribed. There are 
few exceptions, as the papers of Ohtsuka [Ohtl — 3] or Stumpf & Le [Stic]. The 
geometry of the crack can be prescribed in a strong form, like in the case of a plane 
rectangular or elliptic crack which is supposed to remain plane rectangular or elliptic 
during its growth. In a weak form, the geometry of the crack can be prescribed by 
the assumption that the configuration of the body is 2 dimensional and the crack is 
supposed to have only an edge, which is a point. In this case the evolution of the 
crack is conveniently reduced to the evolution of a point. Under these assumptions 
the geometrical nature of the main unknown is obscured. 

A new direction of research in brittle fracture mechanics begins with the article of 
Mumford &; Shah [MS] regarding the problem of image segmentation. This problem, 
which consists in finding the set of edges of a picture and constructing a smoothed 
version of that picture, it turns to be intimately related to the problem of brittle crack 
evolution. In the before mentioned article Mumford and Shah propose the following 
variational approach to the problem of image segmentation: let : C — > [0, 1] 
be the original picture, given as a distribution of grey levels (1 is white and is 
black), let u : CI ^ R he the smoothed picture and K be the set of edges. K 
represents the set where u has jumps, i.e. u G C^{Q\ K,R). The pair formed by 
the smoothed picture u and the set of edges K minimizes then the functional: 

Iiu,K) = [ a \ Vuf dx + [ P \ u-g\'^ dx + jH^{K) . 

The parameter a controls the smoothness of the new picture n, /? controls the 
distance between the smoothed picture and the original one and 7 controls the total 
length of the edges given by this variational method. The authors remark that 
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for P = the functional / might be useful for an energetic treatment of fracture 
mechanics. 

An energetic approach to fracture mechanics is naturally suited to explain brittle 
crack appearance under imposed boundary displacements. The idea is presented in 
the followings. 

The state of a brittle body is described by a pair displacement-crack, (u, K) is 
such a pair if ii' is a crack — seen as a surface — which appears in the body and u is 
a displacement of the broken body under the imposed boundary displacement, i.e. u 
is continuous in the exterior of the surface K and u equals the imposed displacement 
uq on the exterior boundary of the body. 

Let us suppose that the total energy of the body is a Mumford-Shah functional 
of the form: 



The first term of the functional E represents the elastic energy of the body with 
the displacement u. The second term represents the energy consumed to produce 
the crack K in the body, with the boundary displacement uq as parameter. Then 
the crack that appears is supposed to be the second term of the pair (u, K) which 
minimizes the total energy E. 

We shall use this idea in the study of quasi-static brittle crack evolution. For do- 
ing this, we proceed to a time discretization which transforms the problem of crack 
evolution into a sequence of energy minimization problems. Francfort &: Marigo 
[Fma] proceed in the same way in the case of brittle brutal damage evolution. How- 
ever, it is only a belief that when the time step goes to zero, the discretized evolution 
converges to an almost continuous (with respect to time) evolution. We have found 
in the frame of generalized minimizing movements, introduced by De Giorgi [DG], 
stronger mathematical reasons to support this belief. That is why we introduce 
the notion of energy minimizing movement as a particular case of a generalized 
minimizing movement. 

In section 2. the notion of energy minimizing movement is introduced in a form 
useful in the sequel. After preliminaries concerning the statics of a brittle body, the 
Griffith criterion of brittle crack propagation is presented in subsection 3.3., as a 
selection criterion amongst all possible crack evolutions. At the end of this section 
we formulate the problem of quasi-static brittle crack evolution in the form ()3.3.8|) . 

In section 4. we give an energy minimizing movement formulation to this prob- 
lem by using a Mumford-Shah energy functional (definition 4.1.). We are assured 
about the existence of the discretized (or incremental) solution to that problem by 
theorem 4.1. (in a weak form presented in the section of proofs). In subsection 
4.2. we investigate the features of this first model. In this model we have only 
one material constant connected to fracture, namely the constant of Griffith G. We 
find exact solutions and useful estimations in the case of anti-plane displacements 
(theorem 4.2.), which tell us that crack appearance is allowed in this model. A part 
of this results can be found in [Bu2], in connection with fiber-matrix debonding in 
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composites. We prove also a bad feature of the model: the critical stress which lead 
to fracture in an uni-dimensional traction experiment is not a constant of material. 

In section 5. is presented an improved model, based on a preliminary study 
of smooth brittle crack propagation (see also [Bui], [Bu3]). The Griffith criterion 
is reformulated by using proposition 5.1. and the K2 functional (definition 5.3.). 
The K2 functional is a generalized version of the J integral of Rice [R]. After we 
modify the Griffith criterion by the extension of the functional K2, we obtain the 
differential global criterion of brittle crack appearance (DA). A stronger version of 
(DA) is the local crack appearance criterion (LA). In subsection 5.4. the improved 
model is presented. In this model we have two constants of material connected to 
fracture: G and a quantity with dimension of stress named S. The critical stress 
which lead to fracture is deduced from the the elastic constants and E, hence this 
time it is a constant of material. 

Section 6. is devoted to a brief introduction to special functions with bounded 
variation or deformation. Weak versions of theorems 4.1. and 5.2., concerning the 
existence of the discretized (or incremental) solutions of the models presented here, 
are given as consequences of more general results due to De Giorgi k. Ambrosio 
[DGA], Ambrosio [Al— 3], Belletini, Coscia & Dal Maso, [BCDM]. 

In section 7. are given the conclusions regarding the features of the two models 
presented in the paper. We prove a general existence result of the energy minimizing 
movement described in the first model under the assumption of uniformly bounded 
power communicated by the rest of the universe to the body. A comparation is made 
with the model of Ambrosio & Braides [AB], based also on generalized minimizing 
movements, where the crack appearance is forbidden. 

2 General energy minimizing movements 

An energy minimizing movement is a particular case of a generalized minimizing 
movement. The latter notion has been introduced by De Giorgi in [DG], inspired 
by the paper [ATW] of Almgren, Taylor &: Wang. The definition of a generalized 
minimizing movement is (according to Ambrosio [Amb]) the following: 

Definition 2.1. Let S he a topological space and 

F : (1, +oo) xATxS'xS'^RU {+00} 

he a function. For any uq E S, a function u : [0, +00) S is a generalized min- 
imizing movement associated to F with initial data uq, if there exists a diverging 
sequence (sj)jgAr, Sj > 1, and there are functions Ui : N ^ S such that: 
i) Ui{0) = uo; 

a) for any k E N and any i, Ui{k + 1) minimizes the functional 

V F{si,k,v,Ui(k)) 

over S; 
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in) for any t > 0, ttj([sji]) u{t) in S as i ^ +00 . 

The canonical example of (generalized) minimizing movement is given by the 
choice: S = BP, f : —>■ R Lipschitz continuous and and 

F{s,k,u,v) = f{u) + -\u-vf . 

In this case, for any uq G M" there is only one minimizing movement, namely the 
unique solution of the Cauchy problem 

u'{t) = -Vf{u{t)) ,u{0) = uo . 

An energy minimizing movement is a generalized minimizing movement associ- 
ated to a particular function F. It is designed to be a "weak stable" solution of an 
evolution problem of the following type: 

' A{u{t),a{t),t) =0 V t > 

< ^a(t) < L(a(t),n(t)) Vi>0 (2.0.1) 
^ u{0) = Uo ,a(0) = ao 

There arc two unknowns in this problem: u and a. The evolution of these 
unknowns is quasi-static with respect to u. Suppose that we are not in position 
to give a proper law of evolution of a, or that the law of evolution that we have 
gives too many solutions. Assume further that we have instead at our disposal the 
expression of the total energy of the system described by the pair {u,a), name it 
f{u,a), and a set of constraints, not in a differential form, upon the evolution of 
a. We can make then a time discretization with time step S and recursively find 
{uf._^i,af._^-^) from («^,a^), by a minimization process of the total energy / under 
some constraints. A weak stable solution of the previous problem will be a limit of 
sequences (ti^,a|.)fc when the time step 6 converges to 0. 

In the further definition S may be seen as the space of all pairs x = {u,a), 
endowed with a topology. 

Definition 2.2. Let S be a topological space and 

F : {1,+(X)) X N X S X S ^RU{+oo} , 

F{s,k,x,y) = f{s,x,y) + ip{k/s,y) 

be a function, with f: NxSxS^R and : [0, 00) x 5 ^ {0, -foo}. For any 
xq E S, a generalized minimizing movement x : [0, -|-oo) — ^ S associated to F with 
initial datum xq is an energy minimizing movement associated to the energy f with 
the constraints tp and initial datum xq. 

Let us denote by S{\) the following set: 
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S{X) = {yeS: ^{X,y) = 0} . 

Prom definition 2.2. we see that x : [0, +oo) 5 is an energy minimizing 
evolution associated to /, with the constraints ■0 and initial data xq if there exists 
a diverging sequence (si)igjv , Sj > 1, and there are functions Xi : N ^ S such that: 

i) Xi{0) = xq ; 

ii) for any k e N and any i e N , Xi{k + 1) minimizes the functional / over the 
set S{k/si) (in particular Xi{k + 1) belongs to S{k/si)); 

iii) for any t > 0, Xi{[sit]) — > x{t) in S as i ^ +oo . 

3 Notations and preliminaries 

3.1 Notations and constitutive assumptions 

The open bounded set C represents the reference configuration of an elastic 
body and u : $7 ^ is the displacement field of the body with respect to this 
configuration. We shall always suppose, without mentioning further, that the open 
set O and its closure have the same topological boundary. 

The expression of the elastic (or free) energy of the body is: 

w{Vu) dx . 
The first Piola-Kirchoff stress tensor S is 

S(u) = -(Vu) 

and the equilibrium equation of the body in the absence of volumic forces is 

div S(u) = in J) . 

In this paper we consider that the body is linear elastic and homogeneous, i.e. the 
function w{u) has the form: 

w{u) = ^CVu : Vu , 
with the elasticity 4-tensor C having the symmetries: 

^ijkl ^jikl Cklij 

Under these assumptions the stress tensor S becomes the Cauchy stress tensor: 

a = a{u) = CVu = Ce(u) , 



/ 

Jc 
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where e(u) is the symmetric part of Vu, i.e. 

6(u) = ^ (Vu + (Vu)^) . 
We shall suppose moreover that w satisfies the growth conditions: 
VF G IR\ F = F'^ , c I F |2 < w{F) < C | F 1^ , 

where c and C belong to (0, +00). 

In the case of plane displacement the domain f2 C represents a section in the 
cylindrical reference configuration of the body Q x R and u : O — ^ is a plane 
displacement. The displacement with respect to the 3 dimensional configuration of 
the body has the following expression: 

{xi,X2,X3) G U xRt->- {ui{xi,X2),U2{xi,X2),0) . 

In this case we suppose that the body is linear elastic, homogeneous and isotropic. 

In the case of anti-plane displacements the domain 17 C represents a section 
in the cylindrical reference configuration of the body 0, x R too. The anti-plane 
displacement is a function u : ^ R. The 3 dimensional displacement has the 
following form: 

(xi,X2,X3) £ ri X M i-> (0,0,ti(a;i,X2)) G . 

In this case we make the same assumption of isotropic body, therefore the elastic 
energy takes the form: 

fj, I Vu p dx , 
where fi is one of the two Lame's constants. 



If 



3.2 Statics of a fractured elastic body 

For any measurable set C M", | i? |= jC"{B) denotes the Lebesgue measure of 
B and l-i^{B) denotes the k dimensional Hausdorff measure of B. 

By a crack set in the body J7 we mean (according with Ball [Ba]) a topologically 
closed countably rectifiable set, generically denoted by K. 

Given the function /, a point a; G O and an unitary vector (or direction) n G M", 
the approximate limit of / in x with respect to the direction n is denoted by /(x, n) 
and it is defined by the following expression: 

^.^ iBp{x)r\{y : {y-x)-n>o} i f^y) ~ I _ ^ 

P ™+ I Bp{x) n {y : (y - x) • n > 0} I 

Whenever a field of normals at K is chosen, the lateral limits /"•" and /~ of any 
function f : n\K &s:q f+ : K ^ R sind f' : K ^ R, defined by 

f'^ix) = /(x,n(x)) , r{x) = f{x,-n{x)) . 
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This means that and / satisfy the equahties: 

VxEK, hm ^^-(-^^^^ (.-.).n>o}l/(^)-/^(^)l 
I Bp{x) n{y : {y-x)-n>0}\ 

yx^K, hm (.-.).n<o}l/(^)-/-(^)l ^ ^ 

P^o+ I Bp{x) n{y : {y - x) ■ n < 0} \ 

Remark that for any x £ K the object {f^{x), f^{x), n{x)) it is unique to a change 
of signs, i.e. 

(/+(x),r(x),n(x)) ~ (r(x),/+(x),-n(x)) . 

For a given crack set IT in 0, by an admissible displacement with respect to K 
we mean a function u : \ X — > R'^ (where k might be 1,2 or 3) which is and 
posses continuous lateral limits on K. In this section we shall consider the space 
W^''^{Q \ K) n L°°(J7) as the set of weak admissible displacements with respect to 
the crack set K. 

Let n be the dimension of the reference configuration 0. For a given uq G 
{dn, i2")nL°°(90, R"") and for a given rectifiable crack set K, such that W-'^{dn\ 
K) > 0, the following problem has a solution u = u(uo,i^), unique to rigid dis- 
placements o{n\K equals to on 50: 

div ct(u) = in \ i^T 

o-+(u)n = cr-(u)n = on ivT (3.2.1) 
u = uq on \ K 

We use the same notation — u = u{uq,K) — in the anti-plane case, when n = 2, 
k = 1 and the problem ()3.2.1|) becomes 

/i div Vn = in Q\K 

(VM)+n = (Vu)~n = on K (3.2.2) 
u = Uq on d^l \ K . 

The solution u(uo,i^) of the problem (|3.2.H) minimizes the functional 

E{w) = i w{Vw) dx 
Jn 

over the following set of weak admissible displacements with respect to the crack set 
K: 

{v G W^'^{n \ K, i?") n L°°(S7, i?") : V = uo on \ K} . 
By standard arguments it follows that the functional 

V G VF^'2(rj,i?") ^ / cj(u(uo,K)) : Vv dx 
Jn 
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depends only on the trace of v on dQ, hence it give raise to the hnear continuous 
function: 

T(i^) ■.m{dn,H!')nL°^{dn,K')^H-^{dVL,R!') , 

(T(K)uo,v) = / o-(u(uo,K)) : Vv' dx for any v' = v on . (3.2.3) 
Jn 

In the latter definition (•, •) is the duality product of the pair of spaces H^^i^dVL^R^) 
and i^~2 i?"). The function T{K) is called the Dirichlet-to-Neumann map of 
the elastic body $7 with the crack set K. 

Under the assumptions concerning the elastic energy density more precise 
because of the symmetries of the elasticity tensor C, the function T{K) is also 
self-adjoint, i.e. for any u, v we have 

(T(K)u,v) = (T(K)v,u) . 

In the same way can be defined the Dirichlet-to Neumann map associated to the 
problem (HTT^ . 

Remark finally that, under the assumptions considered for w, the elastic energy 
of the body can be expressed using the Dirichlet-to-Neumann map. Indeed, we have: 

/ u;(Vu(uo,K)) dx = i(T(K)uo,uo) . 

3.3 The Griffith criterion of brittle crack propagation 

Let us consider in the elastic body an initial crack set Kq which evolves and 
becomes at the moment t the crack set Kt- We assume that the crack set always 
increase in time, i.e. 

< t < t' , Kt C Kt> . (3.3.4) 

We suppose that the evolution of the body is quasi-static. At the moment 
t the state of the body is characterized by the pair {\i{t),Kt), where u(t) is the 
displacement of the body, admissible with respect to Kt. Let us denote by uo(4) the 
trace of u(i) on dVt. We have then the equality u(t) = u(uo(t), i^t). 

The power given to the body by the rest of the universe at the moment t has 
the following expression: 

P{t) = [ S(u(t))n-uo(0 dx = {T{Kt)uo{t),uo{t)) ■ 
J on 

Let us consider a given curve t {u(t),Kt), such that for any t we have u(t) = 
u(uo{t) , Kt) . For a given t we introduce the following curve of displacements: 

Vt>0, w(r) = u{uo{t + T),Kt) . 
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w(r) represents the displacement the body at the moment i + r in the presence of 
the crack Kt. An easy calculation lead us to the equality: 

^^u;(Vw(r))dx|^^„ = P{t) . (3.3.5) 

Therefore P{t) represents the power consumed at the moment t by the body in order 
to modify its displacement, constrained to follow the path of imposed boundary 
displacements 1 1— > Uo(i), without any modification of the actual crack set Kt. 

The Griffith criterion of brittle crack propagation asserts that during the prop- 
agation of the crack Kt the following inequality is true at any moment t: 

j^y^w{Vu{t)) dx + GH"-H^*)} < Pit) ■ (3.3.6) 

Here G is the constant of Griffith, supposed to be a material constant. 

The relation (|3.3.6p can be written in a different form using the Dirichlet-to- 
Neumann map T{Kt). Let us assume that the crack evolution is smooth in the 
sense that the function t T{Kt) is differentiable, i.e. the Dirichlet-to-Neumann 
map varies smoothly in time. The Griffith criterion takes the following form: 

i(|[T(i^,)]uo(t),uo(t)) + ^{T{Kt)Mt),Mt)) + 

+ ^{T{Kt)uo{t),uo{t)) + Gj^{n^-\Kt)] < {T{Kt)uo{t),Mt)) ■ 

The function T{Kt) is self-adjoint, therefore we obtain the following expression of 
the Griffith criterion: 

i(|[T(KO]uo(t),uo(t)) + G^^{n--HKt)} < . (3.3.7) 

We can see that we have the following equality: 

Pit) - jJ^w{Vu{t)) dx = -i(|[T(i^,)]uo(i),uo(t)) . 

The quantity from the left of the previous equality is usually called the energy release 
rate due only to crack propagation. 

It is obvious that uo(t) plays the role of a time-dependent parameter, since in 
the last inequality uo(*) does not appear. 

The problem of quasi-static brittle propagation of an initial crack in an elastic 
body under a time-dependent imposed displacement uo(t) can be formally put in 
the form H2.U.1() . If we put apart the constraint 1)3. 3. 41) . we have the following 
formulation: 

u{t) - u{uoit),Kt) = Vt>0 
l(|[T(i^O]uo(t),uo(t)) + Gi{n^-HKt)} < yt>0 (3.3.8) 
u(0) = uo , Ko = K . 
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4 The first model 

In the left term of the Griffith criterion appears the time-derivative 

of an energetic functional. Let us consider the following set of admissible pairs 
displacement-crack: 

M = {{u,K) : K is a crack set and u G C^{Tl\K,W) such that 

(u'*',u~,n) exists on . 
The Mumford-Shah energy functional over M has the following expression: 

/:M^MU{+oo} , I{u,K) = [ u;(Vu) dx + . (4.0.1) 

Jn 

4.1 Introduction of the first model 

It is natural to try to give an energy minimizing movement formulation to the 
problem ()3.3.8() by using the functional defined at H4.0.1|) . According to definition 
2.2. and the constraint (|.S.3.4|) . we give the following formulation: 

Definition 4.1. Let us consider the space M endowed with the topology given 
by the convergence: 

We define the functions 

J:MxM^R, J((u,i^),(v,L)) = [ w{Vv) dx + Grr-^{L\K) , 

Jn 



^' : [0, cx)) X M ^ {0, +00} , ^(A, (v,K)) = 




if "v = uo(A) on d^l \ K 
otherwise 



We consider the initial data (uo,i^) G M such that uq = u(uo(0),K). 
For any s > 1 we recursively define (u^, K'^) : N ^ M like this: 

i) (u^i^^)(0) = (uo,K) ; 

ii) for any k £ N {u'^,L^){k + 1) £ M minimizes the functional 
(y,L)£M ^ J(((u^K^)(A;),(v,L)) + ^{{k + 1)/ s, {v, L)) 

over M . In order to verify the constraint S.J^ , K^{k^\) is defined by the formula: 

K'{k^\) = K'{k)\JL\k^\) . 

An energy minimizing movement associated to J with the constraints ^3.3.4}) , ^ 
and initial data (uo,K) is any (u, ii') : [0,+oo) M having the property: there is 
a diverging sequence (sj) such that for any t > 

K'-'){[sit]) -^{u,K){t) as i^ 00 . 
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In the previous definition 1 / s is the step of the discretization of the time variable, 
hence (u*(/c), K^{k)) represents the approximate pair displacement-crack at the time 
k/s. We name any function (u'*, K^) : N ^ M an incremental solution if it verifies 
i) and ii) from the definition 4.1.. 

When Si converges to oo the time step goes to and the incremental solution 
(u*% converges to {u,K){t), for any t > 0. 

A necessary condition for the existence of an energy minimizing movement in- 
troduced in definition 4.1. is that for any given s the incremental solution k G 
N I— > (u*,K*)(fc) exists. The following theorem provides an answer to this exis- 
tence query. In the general case n = 3 this theorem is true, to our knowledge, only in 
a weak form, presented in section 6. In the anti-plane case, however, due to partial 
regularity results for the minimizers of Mumford-Shah functional from [DGCL] , the 
theorem is certainly true. 

Theorem 4.1. Let C M" be a bounded open set with piecewise smooth bound- 
ary, let {uq,K) be a given admissible pair displacement- crack in Q and let 

Uo:N ^ H^dn,R'')nL'^{dn,BJ') 

be a given sequence of imposed displacements such that Uq = u(uo(0),if) ond^l\K. 
Then there exists the sequence (u, K) : N ^ M such that: 

i) u(0) = uo and K{fd) = K; 

ii) for any k & N there is a crack set L{k -\- 1) such that (u(fc -|- l),L{k -\- 
1)) G M, u{k + 1) = uo(A; + 1) on dn \ L{k + 1) and {u{k + l),L{k + 1)) is a 
minimizer of the functional 

{v,L) e M ,v = uo{k-\-l) on dn\L ^ J {{u{k),K{k)),v,L)) . 

The set K(k -|- 1) is given by the formula 

K{k + 1) = K{k) U L{k-\-l) . 



4.2 Features of the first model 

We shall investigate further the behaviour of the model proposed in definition 
4.1. in the particular case of anti-plane displacement. There arc some obvious 
adjustments to be made, is now a bounded domain in and the displacement 
is a scalar function u. The functional J will take the form: 

J{{u,K),{v,L)) = [ fi \ Vv\'^ dx -\- GT-0-{L\K) . 
Jn 
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Let us consider a particular type of imposed displacement on dfl. We split the 
boundary of the body in three parts: 

n = , fi n f2 = , nHri) ■ nHrl) ■ nHTf) > o . 

At any moment t > 0, Tf is force free, i.e. the displacement is not prescribed on 
this part of the boundary. On and the imposed displacement is defined by 
the formula: 

r4^\r \ f oil 
^'^(^^^^^ = [t6 ouf! , 

where 5 is a positive constant with dimension of speed. This displacement is homo- 
geneous with respect to the time variable: 

Vt > , uo{t) = tuo(l) . 

We suppose further that at the moment t = there are no cracks in the body. This 
assumption takes the form K = 0. At f = we have uo{0) = 0, hence the initial 
data are {uo = 0, K = 0). 

Let us consider a time discretization given by the parameter 1/s and the incre- 
mental solution k G N {u^ , K^){k) introduced in definition 4.1. for the initial 
data and the imposed boundary described above. In order to shorten the notations 
we shall omit for the moment the superscript s. 

The incremental solution {u, K) : N ^ M is recursively defined by the following 
two rules: 

i) u{0) = and K{0) = 0; 

ii) for any k G N we seek for the crack set L{k+1) and for the displacement 
u{k + 1) such that {u{k + l),L{k + 1)) G M, u{k + 1) = {k + l)/s uo{l) on 
(ri U Tl) \ L{k + 1) and {u{k + l),L{k + 1)) is a minimizer of the functional 

{v,L)eM ,v = ik + l)/suo{l)on {TiuTl)\L ^ J {{u{k) , K (k)) , (v , L)) . 

The set K{k+ 1) is given by the formula 

K{k + 1) = K{k) U L{k + l) . 

Let us denote by the displacement of the body fi, without cracks, under the 
prescribed displacement on the boundary uq{1). With the use of a notation made 
before, U0 is defined by U0 = u(no(l),0). For any k G N we have (k/s ti0,0) € M 
and k/s uq = k/s Uo{l) on F^ U F^. Therefore, for any k £ N we have 

Ji{u{k),K{k)),{u{k + l),L{k + l))) < J(KA:),i^(A:)),((A: + l)/s «0,0)) . 

The last inequality reads: 

[ fi |VM(fc + l)p dx + Gn^{L{k + l)\K{k)) < 
Jn 



(7) l^^^l' • 

(4.2.2) 
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We can always find a curve in $7 which separates from F^. Moreover, we 
can find such a curve which is a length minimizer in the family of all curves in Q 
separating F^ from F„. Let us denote this curve by S (which exists but it might not 
be unique). The domain Q has the following decomposition with respect to S: 

n = f)^ u 0^ , F^ c , F^ c , n = , 
Wnf^ = s . 

Let us define the following displacement: 



usix) 



x£Q^ 
6 X G ^2 



It is easy to see that for any k £ N we have {k/s us, S) £ M and k/ s ug = kj s uo(l) 
on (F^ U F^) \ S. Therefore we obtain the following inequality: 

/ H |Vu(A; + l)p dx + G'H}{L{k^\)\K{k)) < Gn\S\K{k)) . (4.2.3) 
Jn 

From (|4.2.3() we derive the following conclusion: for large time k/s the crack 
set K{k) is not void. Indeed, suppose that the function k £ N i-^ {k/s U0,0) 
is an incremental solution constructed by the rules i) and ii) above. Then for any 
k £ N the inequality 1)4.2. 2() becomes an equality and the inequality 1)4. 2. 3|) takes 
the following form: 

y /i I Vn0 |2 dx < Gn\S\K{k)) , (4.2.4) 

which lead to contradiction. Therefore this model can predict crack appearance. 
The critical step k, after which a crack appears in the body (as the incremental 
solution predict), is the greatest natural with the property 1)4. 2. 4|) . 

The following theorem contains stronger informations regarding the minimizers 
of the Mumford-Shah functional in our particular case. 

Theorem 4.2. Let $7 C be a bounded open set with piecewise smooth bound- 
ary dil. and let n be the field of outward normals over the boundary. Let us suppose 
that the boundary of Q has the following decomposition: 

FlnF^ = 0,finf2=0, nHri) • ^^(f^) • n'Tj) > o . 

Let us consider the functional 

L{v,K) = I [ \ Vv\ dx + Gn^{K) , 
2 Jn 
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defined over the set 

{{v,K) : V e C^{Tl\K,R)} . 
Let, for any D G R, u{D) -.Vt^Rhe the solution of the problem: 

div = in ^} 
Vvn = onTf 
V = on 
^ V = D on . 

We suppose that exist strictly positive numbers c and C such that for any x G U 

C >| Vti(l)n \{x)>c . 

There exist then two numbers m < M , which depends only on r2,F^,F^ and Vf, 
such that: 

i) if < m then (ti(D),0) is the only minimizer of the functional I over the 

set 

M{D) = {{v,K) : V eC^(Sl\K,R) , V = u{D) onVi^JVl] , 

ii) if > M then any minimizer of the functional I over the set M{D) has the 
form (uk,K), with \ Vuk |= almost everywhere in Q and K geodesic in (i.e. 
length minimizer) separating F^ from F^. 

Moreover, if c = C then M = m, hence if VM(l)n is piecewise constant on 
F^ U F^ then we have only two kinds of minimizing crack sets. 

The theorem assures us that for small time k/s the body remains sane and for 
large time k/s a crack with a particular shape appears in the body. Precisely, 
for small k/s we have {u{k),K{k)) = {k/s ii0,0) and for large k/s we have 
{u{k),K{k)) = {k/s us,S). The theorem help us to find particular cases when 
the passage from the first type of minimizer to the second one is brutal. Indeed, 
consider that is a rectangle (0, a) x (0,L), 

Tl = (0,a)x{0} , ri = (0,a)x{L} 

and Fj is the remaining part of the boundary. Let us consider, for simplicity, that 
5 = 1. With the notations from the theorem 4.2. we have c = C therefore we 
have only two kinds of pairs displacement-crack which compete. We use ()4.2.4() in 
order to find the critical time k/s when the incremental solution {u,K) : N ^ M 
switches from {k/s n0,0) to {k/s us,S), where S is, for example, (0,a) x {L/2}. We 
find that the critical A;/s is determined by the double inequality: 

kV ^ GL ^ fk + l^^ 



s J ji \ s 

We are lead to the definition of the critical moment tc, given by the formula 

,2 _ GL 
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tc is proportional with the square root of L. The anti-plane stress existing in the 
sane body at the moment tc has the following expression: 




We can see that this stress depends on L, hence on the geometry of the body. 
Because G is supposed to be a material constant we obtain the following conclusion: 
the model described above is not compatible with any model of crack appearance based 
on a critical stress as material constant. 

5 The improved model 

We have seen that the first model allows crack appearance but it is not compatible 
with any critical fracture stress based model. Our purpose is to improve the first 
model in order to allow the existence of a critical stress which damages a structure. 
We shall find a way to make this improvement by studying first how smooth brittle 
propagation of cracks can be described with the Mumford-Shah energetic functional. 

5.1 Smooth brittle crack propagation 

There are two steps in order to define the notion of smooth brittle crack prop- 
agation. The first step consists in smoothness demands on the initial crack set K. 
We shall suppose that K is endowed with the structure of manifold with boundary. 
The boundary of K, denoted by dK, represents the edge of the crack. The second 
step consists in smoothness demands on the evolution t ^ Kf of the crack. We shall 
restrict our attention only to evolutions of the initial crack K obtained by smooth 
deformations of K. The initial crack may be as complex as we wish, because the 
structure of manifold with boundary allows that, but this complexity remains the 
same during the propagation of the crack. 

We shall work with deformations of the initial crack set K by endomorphisms of 
J7. Let us consider the following set of diffeomorphisms: 

v = c°°{n, n) n w''^{n, i?") : (p-^ e c^in, o) and supp (</> - in) c n} . 

(5.1.1) 

We have denoted by n the dimension of the space where lies. The introduction 

of the Sobolev space W^''^{fl, R^) has been made for mathematical reasons ( to 
be found for example in Ebin & Marsden [EbM]) and the number s is chosen to 
be greater than § + 2. In the paper [Bu3] a rigorous mathematical description of 
smooth brittle crack propagation can be found. We mention that the number s 
controls the variation of the smoothness of the deformed crack set (f){k) from the 
smoothness of the initial crack set K. 

The condition supp (cp — la) C O, means that near the boundary of O equals 
the identity map. 
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Definition 5.1. A smooth fracture curve is a function 



t G [0, T] ^ (l)teV' , 

which has the following properties: 

i) (Ao = In , 

ii) the map t € [0,T] (/)t G is continuous with respect to the 
topology induced by the norm max (|| • || • llvK'''^); for every t E [0,T] cpt exists 
and 

rjt = cj)t4T^ e W''\Q,R'')nL'^{n,R'') , 

Hi) for any t < t' we have (ptiK) C (f)t'{K) . 
We have used the notation f.g for the composition of the function f with g. 
A crack evolution curve is associated to the smooth fracture curve t ^ (j)t and 
initial crack K by the formula: 

Kt = UK) . 



There are infinitely many smootli crack propagation curves t ^ (j)t witli tlie same 
associated crack evolution curve t ^ Kt- 

Under smoothness assumptions on the initial crack set K, for a smooth crack 
propagation curve t ^ (f)t the condition that the crack grows implies that for any 
t > we have: 

( Vt-n = on (j)t{K) 

I lMK)diVsVtdn--' >0 . ^^■'■') 

The integral from the last inequality equals the variation of the area of the crack set 
(see Allard [All]). The operator divs is the tangential derivative with respect to the 
surface (ptiK) and it has the following form: 

divst] = div r] — n ■ (Vry) n , 

where n is the normal to the surface (j)t{K). 



5.2 K2 functional and the Griffith criterion 

We want now to reformulate the Griffith criterion of brittle crack propagation 
()3.3.6|) in terms of smooth crack propagation curves. Our assumptions on the evo- 
lution of the body are the following: 

Al) the evolution of the linear elastic body is quasi-static, 

A2) a smooth curve of imposed displacements 1 1— > uo(t) is given on 

the boundary 

The assumption A2) can be modified by the replacement of dO. with a fixed part 
of the boundary r„; on the remaining part we suppose that the body is force 
free. 
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At any moment t the state of the body is described by the admissible pair 
displacement-crack {u{t),Kt). The assumption Al) implies that the displacement 
u{t) is determined by the knowledge of Kt and boundary condition uo(t). With 
a notation used several times before, we have u(i) = u{uq (t), (j)t{K)). We shall 
change for our purposes this notation by writing: 

u(t) = u(uo(t),0f) . 

Definition 5.2. Let t e [0,r] ^ u°(t) G C(5J7,i?'') be a curve of imposed 
displacements on the exterior boundary of the body. A balanced fracture curve is any 
function 

t e [0, T] ^ (u*, (l)t) G W^'^in \ K, i?") X V , 

satisfying the following items: 

i) for any t G [0, T] we have 

Ut-'Ar^ = = u(uo(t),(/)t) , 

ii) t ^ (j)t is a smooth crack propagation curve. 

For given curve of boundary displacement and initial crack set for any smooth 
crack propagation curve there is only one associated balanced fracture curve. The 
Griffith criterion of brittle crack propagation will act as a selection criterion amongst 
all smooth crack propagation curves. 

We have seen that an equivalent form of the Griffith criterion is (|3.,3.7|) . With 
the change of notation T{(j)t) = '^{4't{K)), we say that a smooth crack propagation 
curve is compatible with the Griffith criterion if for any t we have 

i(|[T(<At)]uo(t),uo(t)> + Gj^{H^-\UK))] < . (5.2.3) 

The first term from (|5.2.3() represents the variation of the elastic energy of the 
body calculated for the following variation of the displacement: 

T ^ U(t + T) = u(uo(t),0t+r) • 

The dependence of u(t + T) with respect to (j){t + T) is implicit. An explicit variation 
of the displacement would be preferable, like this one: 

v(r) = M{t).<t)t4tlr ■ 

We are lead, by a change of variables, to the following equality: 

2 ( / w{Vw{t)) dx ) = / {[CVu(uo(t),0t) : Vu(uo(t), dii^ m - 
dr \Jn ) Jn 

(5.2.4) 
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-2 [CVu(uo(t), [Vu(uo(t), <t}t)]ik [Vr?4j dx . 

Proposition 5.1. Lett ^ (pt be a smooth crack propagation curve, rjt = <f>t-<Pi^^ 
and uq £ Hi (^d^l)r\L°°{dQ). Let us define, for fixed t, v(r) = u{t).(pt-4't+T ■ Then 
the following inequality is true : 

(|to,)]uo,uo) < 2 A |^^„(Vv(T))dx|^^„ . (5.2.5) 

This proposition, togetlier witli tlie equality (|5.2.4|) . allows us to introduce a 
generalization of the J integral. We use the notation r] £ W^' {n,R'') for r] £ 
WQ''^{n, i?") with nuh trace on dn. 

Definition 5.3. The generalized J integral is the following functional 
K2 : P"x|Wo'^(l^,ii") nL°°(J^,i?")}x|if5(5J^) n L°°(917, ii")} mathbbR , 

K2{<p,rj,uo) = ^[CVu(uo,</>) : Vu(uo,<A)]dit;r?+ (5.2.6) 

+ [CVu(uo, </.)],, [Vu(uo, <f))]ii, [Vv]kj} dx . 

In order to explain why the functional K2 is the generalization J integral, we 
begin by a temporary introduction of easier notations: 

cr = CVu(uo,(/)) ,u = u(uo,(/>) , t« = ^CVu(uo,(/') : Vu(uo,(/>) . 

We define the tubular neighbourhood, of radius r, of the edge d(j){K) of the crack 
set K: 

Br = Br{d(t){K)) = U^^9^^K)B{X, r) . 

The field of normals over d Br{d(f){K)) will be denoted by v, without specifying the 
parameter r. 

If u belongs to then we have 

w ^,1 - c^iiUj^fc rfk^j = [w rji - aiju^k Vk] j - crkm^k,mi Vi + crii^l,ki Vk + crii,i^l,k Vk ■ 

(5.2.7) 

According to the assumption Al), the divergence of the stress field a equals 0. 
We integrate the equality (|5.2.7|) over 0,\Br and we obtain: 

K2((j),rj,uo) = - lim / [wT]i - cj/jU/^fc %] dx . 

'■-^0 Jn\Br 
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By a flux-divergence formula, we are lead to the following expression of K2: 

i^2((;^>, ?7, uo) = lim <^ / {-wq ■ v + anui^k Vk I'i} > + [w]r] ■ ndW~^ . 

[J dBr{d(l>{K)) J J4,{K) 

The functional K2 is interesting in the case when • n = on (j){K), as ()5.1.2() 
suggests. In this case we have: 

K2{^,r],uo) = lim J / - {wr]- ly - aim^k m ^i] dTT^^ \ . (5.2.8) 

'■^0 yJdBr{d^{K)) J 

Let us consider that we are in the case of plane displacements (hence n = 2) and 
that the crack set (p{K) lies on the Oxi axis. If we take rj equal to (1,0) in a 
neighbourhood of the edge of the crack then we have: 

/i:2((/), r?, uo) = lim <^ / - {wvi - akiUk,! fi} dW"M . 

''^0 yJdBr{d^{K)) J 

We recognize in the right term of the equality above the expression of the classical 
J integral. 

We propose the following selection criterion for smooth crack propagation curves: 

A smooth crack propagation curve t ^ (pt satisfies the generalized 
Griffith criterion if at any moment t > we have • n = on (ptiK) and 

K2i^t,Vt,Mt)) > Gjn^'-HMK)) . (5.2.9) 



Ohtsuka [Ohtl — 4] proves that under stronger smoothness assumptions on K 
and on the curve t ^ Kt always exists a smooth crack propagation curve t ^ (j)t 
such that for any t and with our notations we have: 

(^[T((/>t)]uo(t),uo(t)) + K2{4>uVuMt)) = . 

For this reason we consider that (|5.2.9|) is not too strong with respect to the classical 
Griffith criterion. 



5.3 Extension of K2 and admissible cracks 

We want to extend the Griffith criterion of brittle fracture propag ation (l5T9l) 
in order to allow crack appearance. The leading idea is to consider crack evolution 
curves t^Kt which are limits of crack evolution curves of the form t ^ (f)t[K). 

Let t I— > 0^ be a smooth crack propagation curve. At any moment t the vector 
field r]t = (pt-'PT represents the propagation speed of the edge d(j)t{K) of the crack 
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(pt{K). Precisely the restriction of 7]f to d(ptiK) represents the distribution of speed 
of propagation of the points belonging to this n — 2 surface. The appearance of a new 
crack at the moment t is seen as a limit of processes of smooth crack propagation, 
when the distribution of speed rjt develops jumps. 

We shall consider therefore a sequence {r]h)h in the space W^''^{i^, i?")nL°°(r2, i?"). 
For each h we define the following flow t ^ (p^^: 

(/.[f = In , (A^^ = In + Tr]h . 
For small times r we have (p!}: G and for any h we see that 

Let us suppose that r]h converges almost everywhere to rj. Then for any r (j)'^ con- 
verges almost everywhere to (j)r = In + tij. We make the following supplementary 
assumptions: 

51) for small r (pr is almost everywhere injective, 

52) \ has finite Hausdorff n — 1 measure, 

S2) let us denote by Sri the set where rj has no approximate limit 
(the complementary of the Lebesgue set of rj, see for this the section of proofs) and 
by I Df I the total variation measure associated to the distributional derivative of 
the function /; then we have 

I D'l] \ {n\Sn) = . 

The assumption SI) assures us that (jy-r almost everywhere maps different points 
from in different places in 0t-(J7); S2) prevents the case where in the limit appear 
holes in the configuration 4>ri'^) and S3) is a more sophisticated condition which 
says that no strange Cantor sets appear in cpri^)- 

Under these assumptions it is easy to prove that on 5^ the jump of rj satisfies 
the relation: 

[rj] ■ n = . 

Indeed, suppose that in a neighbourhood of a; G S*,, we have [r/] • n < 0. Then the 
assumption S2) is contradicted because a solid neighbourhood of x is transformed 
by (pr in a neighbourhood with a hole, when r > 0; if we have • n > then SI) 
is contradicted because even if local injectivity is respected, the global injectivity in 
the form SI) is not. 

We shall consider therefore pairs (r/, N) where is a topologically closed count- 
ably rectifiable set, rj has tangential jumps on the surface and satisfies the smooth- 
ness assumption rj S Wq''^{Q \ N, i2") n L°°(J7, i?"). The subscript in the notation 
Wq'^{^ \ N, i?") means that = on dQ in the sense of traces. 

Let us suppose that there is no initial crack in the body: K = 9. We perform 
the same calculation for K2(ln,ry, uq) as we did after definition 5.3. and we obtain 
from (|5.2.7)) the expression: 

K2{ln,r],uo) = - [w rji - aijUi^k rjk] ■ . 
Jn 
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Let us suppose that is a surface with boundary and let be a tubular neigh- 
bourhood of dN, of radius r. Because of the assumption [r/] • n = on N, we 
obtain: 

Jn 

It is natural to try to modify the Griffith criterion (|5.2.9|) in order to have a control 
on the integral from above. We propose the following differential criterion of crack 
appearance (DA), which make a selection amongst all crack sets which can appear 
in the body. The constant S, with the dimension of a stress, which appears in this 
criterion, is postulated to be a constant of material. 

(DA). Let us consider the elastic body 0, and the imposed boundary displacement 
uq. a crack set N can appear in the body if there exists a vector field rj G Wq'^{Q \ 
N, R") n L°°(0, R^), such that [ry] • n = Q on N and 

[ aii(uo,0)niUi,fc(uo,0)[%] dW"-i > ||77||lco E ^"-^(A^) . (5.3.10) 
Jn 

We give also a criterion of local crack appearance (LA). This criterion tells us if 
in the point x G O a crack with normal n can appear. 

(LA). Let us consider the elastic body 0, and the imposed boundary displacement 
Uq. In the point x £ Q a crack with normal n can appear if 

sup {aii{uo, 0)niUi,fc(uo, 0)i/fe : g , | zy |= 1 , n • i/ = 0} > S . (5.3.11) 

From ()5.3.1UI) and H5.3.11|) we see that if N is smooth enough and if for any 
X £ N, the criterion (LA) is satisfied for the pair (x, n{x)), where n(x) is the normal 
to N at X, then satisfies the global criterion (DA). 

Let us suppose that the body is a cylinder to x [0, L] and uq imposed on the 
top and bottom of this cylinder such that 

u(uo,0)(xi,a;2,X3) = (0,0,0x3) ,a>0 . 

The stress o"(uo,0) has the form: 

/ \ 
a(uo,0) = , 
V a / 

where a = Ea. If we denote by a the angle between n and the Ox^ axis, we have: 
sup {crii(uo,0)niUi,fe(uo,0)i^i : 1/ G M" , | |= 1} = t^o"^ sin{2a) . 
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The maximum value of this expression is attained for a = tt/4. Therefore in the 
experience of uniaxial traction the (LA) criterion affirms that a crack can appear if 



—a^ > S , (5.3.12) 

and if we have equality in the relation above then the normal of the crack predicted 
by (LA) makes the angle 7r/4 with the axis of the cylinder. The relation (|5.3.12() 
gives us the value of the critical stress for uniaxial traction (which is a constant of 
material this time): 

o-cr = V2ET. . (5.3.13) 

5.4 The improved model 

In this section we propose an improved energy minimizing movement formulation 
to the problem of brittle fracture evolution (|3.3.8|) . The model is based on two 
constants of material connected to fracture, namely G and S previously introduced. 
In this formulation the critical stress which lead to fracture in a traction experiment, 
defined by (|5.3.13|) . is a constant of material. 

We denote by Sn the set of all i' G M" with | |= 1. Let us define the following 
function: 

/oo : mathhhBJ''''' x M"^'^ x 5„ ^ M U {+00} , 



/oo(CT,F,n) 



G if sup {cTiiWiYi^kV}^ : u € M" , |i^|=l, n-z^ = 0} > S 
+00 otherwise . 

(5.4.14) 

The physical dimension of /oo is the same as the one of G. 

With the use of the function /oo the criterion of brittle crack appearance (LA), 
takes the form: given the imposed boundary displacement Uq, a; G O and n G Sn, a 
crack of normal n can pass by x if: 

/oo(cT(uo,0),Vu(u°,0),n) < +00 . 

The idea of the improved model is to consider only pairs displacement-crack (v, L) S 
M admissible with respect to (LA). 

Definition 5.4. Let us consider the space M of admissible pairs displacement- 
crack endowed with the topology given by the convergence: 

iu,,K,) ^ (u,K) ^f I "^u^'j^Jk) ^ . 
We define the functions 

Joo : M X M ^ M , 
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Joo((u,i^),(v,L)) = f w{Vw) dx + I f^{a{u),Vu,ii) dH^-\x) , 

Jn J L\K 

^ : [0, oo) X M ^ {0, +00} , ^(A, (u, K)) = | ^^^^^J'^ [ 

We consider the initial data {vlq^K) G M such that uq = u(uo(0),i^). 
For any s >1 we recursively define (u^, K'^) : N ^ M like this: 

i) K,i^^)(0) = (uo,K) ; 

ii) for any k & N {u'^,L^){k + 1) G M minimizes the functional 
{^,L)eM ^ J^{{{u',K'){k),{^,L)) + ^((A: + l)/s,(v,L)) 

over M . K^{k + 1) is defined by the formula: 

K'{k + l) = K'{k)UL'{k + l) . 

An energy minimizing movement associated to J^o with the constraints \cl.cl.4\ ), 
^ and initial data (uo,i^) is any (u,K) : [0, +00) M having the property: there 
is a diverging sequence (sj) such that for any t > 

K'^){[sit]) -^{u,K){t) as i^ 00 . 



We are interested if for fixed s an incremental solution exists. There is no result 
to our knowledge that assures the existence of a minimizer of the functional 

(v,L)GAf^ [ w{Vv) dx + [ /oo((T(u),Vu,n) d?t:"-^(x) 
Jn J L\K 

in our case. That is why we prefer to modify the function f^. This function imposes 
a cost equal to +00 to the pairs displacement-crack which are not compatible with 
the (LA) criterion. We shall demand a finite but great cost instead, hoping that non 
admissible pairs will not enter in competition with admissible ones. 

Let us consider a number C > G (with the same physical dimension as G) and 
a function 

fc : mathhbR^'''' x R"^" x M" ^ R , 

with the following properties: 

i) fc is positively 1-homogeneous with respect to the third variable, 

ii) for any n € M" such that | n |= 1, if 

sup {uiiHiYi^k^k :i^GM",|i/|=l,n-i/ = 0} > S 

then/c(cT,F,n) = G, 

iii) for any n G M" such that | n |= 1, if 

sup {uiiTdiYi^k^k : V , \ V \= I , w ■ V = < S 
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then /c((7,F,n) > G, 

iv) for any {a, F, n) G M"^" x R"X" x Sn we have F, n) < C . 
As for theorem 4.1., the following statement has been proven to be true only in 
a weak sense, described in the next section. 

Theorem 5.2. Let C be a bounded open set with piecewise smooth bound- 
ary, let (uo, K) be a given admissible pair displacement- crack in and let 

UQ-.N^ H^dn,K") n L'^{dn,R'') 

be a given sequence of imposed displacements such that Uq = u{uq{0),K) ondVL\K. 
Let us consider the functional Jc 

Jc : M X M , 
Jc{{n,K),{v,L)) = [ w{Vv) dx + / /c(a(u),Vu,n) , 

Jo, JL\K 

where fc is chosen to satisfy the assumptions i) — iv) from above. Then there exists 
the sequence (u, K) : N ^ M such that: 
i) u(0) = Uo and K(0) = K; 

a) for any k E N there is a crack set L{k + 1) such that (u(A; + l),L{k + 
1)) e M, u(fc + 1) = uo(A; + 1) on dfl \ L{k + 1) and {u{k + l),L{k + 1)) is a 
minimizer of the functional 

{v,L) e M ,v = uo{k-\-l) ondn\L ^ Jc {{u{k),K{k)),v,L)) . 

The set K{k + 1) is given by the formula 

K{k + 1) = K{k) U L{k-\-l) . 

6 Proofs 

6.1 Weak versions of theorems 4.1. and 5.2. 

This section is dedicated to a brief voyage trough the spaces SBV and SBD. 
The weak forms of theorems 4.1. and 4.2. are direct applications of results listed 
below. 

The space SBV(J7, i?") of special functions with bounded variation was intro- 
duced by De Giorgi and Ambrosio in the study of a class of free discontinuity prob- 
lems ([DGA], [Al], [A2]). For any function u e L^{n,R'') let us denote by Du 
the distributional derivative of u seen as a vector measure. The variation of Du 
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is a scalar measure defined fike this: for any Borel measurable subset B of 0, the 
variation of Du over B is 

{oo 
I Du{Ai) I : U^^i AiCB , AinAj = il} j 
i=i 

A function u has bounded variation if the total variation of Du is finite. We send the 
reader to the book of Evans &i Gariepy [EG] for basic properties of such functions. 
The space SBV(J7, i?") is defined as follows: 

SBV{n, i?") = {u G L\n, R'') : \ Du \ (Q) < +og , \ D'u \ (O \ Su) = O} . 

The Lebesgue set of u is the set of points where u has approximate limit. The 
complementary set is a negligible set denoted by Su- If u is a special function 
with bounded variation then Su is also a (i.e. countably) rectifiable. 

From the Calderon & Zygmund [CZ] decomposition theorem wc obtain the fol- 
lowing expression of -Du, the distributional derivative of u G SBV(f2,i?"), seen as 
a measure: 

Du = Vu(x) dx + [u] n dH''^^ . 

We shall use further the notation /x ^ A if the measure jj, is absolutely continuous 
with respect to the measure A. 

Let us define the following Sobolev space associated to the crack set K (see 
[ABF]): 

W^^ = |u e SBV(J2,i?") : J I Vu p dx + y" [u]^ d'H'''^ < +00 , | D'u |< 
It has been proved in [DGCL] the following equality: 

w^''^{n\K,W)nL°°{n,EJ') = w^'^{n,R'')nL°°{n,EJ') . (6.1. i) 

Therefore if u = u{uq,K) and uq G L°°(t}0,i?") then u is a special function with 
bounded variation. 

A similar description can be made for the space of special functions with bounded 
deformation SBD(r2) can be found in Ambrosio, Coscia & Dal Maso [ACDM]. For 
any function u G L^{Q,,K^) we denote by Eu the symmetric part of the distribu- 
tional derivative of u, seen as a vector measure. We denote also by Ju the subset 
of Q where u has different approximate limits with respect to a point-dcpcndcnt 
direction. The difference between Su and Ju is subtle. Let us quote only the fact 
that for a function u G SBV(0, R"^) the difference of these sets is 'H"~^-negligible. 

The definition of SBD(fi) is the following: 

SBD(0, i?") = {u G L\n, -R") : \ Eu \ (n) < +00 , \ E'u \ (O \ Ju) = O} . 

If u is a special function with bounded deformation then Ju is countably rectifi- 
able. We have a decomposition theorem for SBD functions, similar to Calderon &; 
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Zygmund result applied for SBV functions. The decomposition theorem is due to 
Belletini, Coscia &: Dal Maso [BCDM] and asserts that 

Eu = e(u)(x) dx + [u]0ndWp^ . 

Here means the symmetric part of tensor product and e(u) is the approximate 
symmetric gradient, hence the approximate limit of the symmetric part of the gra- 
dient of u. 

In order to give weak versions of theorems 4.1. and 5.2. let us weaken first 
the space M of pairs displacement-crack. We introduce the new set of weak pairs 
displacement-crack A4: 

M = {{m,K) : K is cj-rectifiable, u e SBD(S7) and | E'm | (1^ \ K) = 0} . 

(6.1.2) 

Given (u, K) G the set K is countably rectifiable but it is not necessarily closed; 
we impose also weaker conditions on the regularity of the displacement u. A direct 
consequence of (|6TT|l is that any pair displacement-crack (u, K) such that u G 
L°°(0,i2") belongs to the set M. 

Let us define the functional J , the weak version of the functional J introduced 
at definition 4.1.: 

J : MxM^^, J{{u,K),{v,L)) = [ u;(e(v)) dx + Gn''-\L\K) . 

Jq 

Before we introduce the correspondent of the function ^ from the same definition, 
let us explain what we mean by u = uq on the boundary of 0,. We consider, for 
simplicity, that uq : dO, R*^ is a continuous and therefore bounded function. Then, 
for any u G SBD(i7), u = uq if the approximate limit of u equals uq in any point 
of do, where the first exists, i.e.: 

\/x £ do,, if 3v{x) such that 
;• /Bp(x)nn I "(^) - ^(^) I dj/ 

limp^O+ — \ o ( \ r.r,\ =0 thenv(x) = Uo(x) . 

Let us consider a curve of imposed displacements A ^ uo(A) G C{dQ, R"). The 
function Q, introduced instead of is defined as follows: 



e(A,(u,i^)) 



e : [0, +00) X ^ {0, +00} , 

if u = Uo and W^-^ {K \ J u) = 
+00 otherwise . 



Definition 6.1. (weak version of definition 4.1.) Let us consider the space M 
endowed with the topology given by the convergence: 



{uh,Kh) ^ iu,K) if 



Uh L"^ ^ u , 
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Let us consider also the function J , the curve of imposed displacements t ^ uo(t) 
with the associated function Q and the initial data {uq,K) S M such that uq = 
u(uo(0),i^). 

For any s > 1 we recursively define (u^, K^) : N ^ M like this: 

i) (u^i^^)(0) = (uo,K) ; 

ii) for any k £ N {u'^,L^){k + 1) £ M minimizes the functional 
{v,L)€M ^ Jii{u',K^)ik),iv,L)) + e((fc + l)/5,(v,L)) 

over M. In order to verify the constraint ^3. 3.J^ , K^{k^\) is defined by the formula: 

K'{k + l) = K'{k)U3usik+i) . (6.1.3) 

An energy minimizing movement associated to J with the constraints ^3. 3.4^ , Q 
and initial data {uo,K) is any {u, K) : [0, +oo) — > Ai having the property: there is 
a diverging sequence (sj) such that for any t > 

{u'\ K''){[sit]) -^{u,K){t) asi^oo . 



Let us remark that the disappearance of the set L'^^k + 1) from the crack-growth 
condition (|6.1.3|) is only apparent, because if {u'^ , L^){k + 1) minimizes the functional 

{v,L)eM ^ Jii{u',K-'){k),iv,L)) + e((A; + l)/s,(v,L)) 

then e((A; + l)/s, (u^L")(A; + l)) = 0, hence 

n^~\K \Ju) = ■ 

In [ACDM] has been proven that functionals like are inferior semi-continuous 
and coercive, hence on closed subspaces V of SBD(r2) the functional 

vGV^ J((K,i^^)(fc),(v,Jv)) 

has a minimizer. Such a closed subspace of SBD(r2) is the space of all weak dis- 
placements V with V = Uq, where uq is a given boundary displacement. Therefore 
the following theorem is true by a trivial induction: 

Theorem 4.1. (weak version) Let C M" be a bounded open set with piecewise 
smooth boundary, let {uq,K) be a given admissible pair displacement- crack in and 
let 

be a given sequence of imposed displacements such that vlq = u(uo(0),i^) ondVl\K. 
Then there exists the sequence (u, K) : N ^ M such that: 
i) u(0) = uo and K(0) = K; 
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a) for any k £ N there is a countahly rectifiable set L(k + 1) such that 
{u{k + l),L{k + 1)) G M, vL{k + 1) = uo(A; + 1) on dn and {u{k + 1), L{k + 1)) is 
a minimizer of the functional 

(v,L) G 7W ,v = uo(/e + 1) on 017 ^ J ((u(A;), i^(A;)), v, L)) . 

The set K{k + 1) is given by the formula 

K{k + 1) = K{k) U Ju(fc+i) • 



The weak version of theorem 5.2. is obtained in the same way. We start by 
relaxing the functional Jc to the functional J^c- 

Jc{{u,K),{v,L)) = [ w{e{v))dx + / /cKu), Vu, n) dW^-^x) . 

Jn Jl\k 

Here cr(u) = Ce(u) and e(u) is the approximate symmetric gradient of u. 

We have the following definition of an energy minimizing movement associated 
to with the usual constraints: 

Definition 6.2. (weak version of definition 5.4. adapted for J'c) Let us consider 
the space M endowed with the topology given by the convergence: 



{uh,Kh) ^ {u,K) if 



\ih L"^ ^ u , 



Let us consider also the function Jc, the curve of imposed displacements 1 1— > uo(t) 
with the associated function Q and the initial data {uq,K) G M such that uq = 
u(uo(0),i^). 

For any s > 1 we recursively define (u*, K'^) : N ^ M like this: 

i) K,i^^)(0) = (uo,K) ; 

ii) for any k £ N (u*,L^)(A; + 1) £ Ai minimizes the functional 
{v,L)eM ^ Jcii{u',K'){k),{v,L)) + e((fc + l)/s,(v,L)) 

overAi. In order to verify the constraint icl.cl.4\ ), K'^^k + l) is defined by the formula: 

K^k + l) = K^k) UJusik+i) ■ 

An energy minimizing movement associated to Jc with the constraints \cl.cl.4\ ), 
Q and initial data {vlq, K) is any {u, K) : [0,+oo) M. having the property: there 
is a diverging sequence (sj) such that for any t > 

{vi'\K'%[sit\) -^{u,K){t) asi^oo . 
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The assumptions i) — iv) on fc from the previous section ahow us to apply the 
main existence result from [ACDM] to the functional jTc- We have therefore the 
following weak version of the theorem 5.2.: 

Theorem 5.2. (weak version) Let Q C M" be a bounded open set with piecewise 
smooth boundary, let {uo,K) be a given admissible pair displacement- crack in and 
let 



be a given sequence of imposed displacements such that Uq = u(uo(0),K) ondQ\K. 
Then there exists the sequence (u, K) : N ^ M. such that: 



i) u(0) = uo and K{0) = K; 

a) for any k £ N there is a countably rectifiable set L{k + 1) such that 
(uik + 1), L(/c + 1)) G M, u{k + 1) = uo(A; + 1) on dn and (u(/c + 1), L{k + 1)) is 
a minimizer of the functional 

{v,L)eM ,v = uo{k + l) ondn ^ Jc ((u(A:), v, L)) . 

The set K(k + 1) is given by the formula 



In the anti-plane case we have to consider the space SBV(i7,i?) instead of 
SBD(ri, i?"). In this case the partial regularity results of De Giorgi, Carriero k. 
Leaci [DGCL] and Ambrosio [A...] tell us that the classical Mumford-Shah func- 
tional has minimizers in the set of pairs displacement-crack M. Both theorems are 
therefore true in the strong form, in the anti-plane case. 

6.2 Theorem 4.2. 

For any D > let uq{D) be the following boundary displacement: 



For any crack set K we denote by uk the displacement uk = u{uq{1),K), i.e. the 
solution (or one of the solutions) of the problem: 



Uo : ^ C{dQ, i?") 



K{k + 1) = K{k) U Ju(fc+i) . 




divVu = in n\K 



Vun = on TfUK 
u = on 
^ u = l on Tl . 
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More general, we shall use the notation uk{D) = u{uq{D)^K). It is obvious that 
uk{D) = D UK- In the body of Theorem 4.2. we have introduced the Mumford- 
Shah functional: 

I{v,K) = I [ \Vv\ dx + Gn\K) . 

The displacement uk has the minimum property 

IiuK,K) < I{v,K) , W{v,K)eM , V = uo(l)onriur2 . 

That is why it is reasonable to redefine the functional J as a functional depending 
only on the crack set K: 

i{K) = I{UK,K) . 
With this notation we have, for any D > 0, the inequality: 

I [ \ VuK I dx + Gn\K) < I{v,K) , 'i{v,K) £M ,v = uo{D) on T^ur^ . 
We make the notation: 

i{K,D) = I [ \ Vuk\ dx + Gn\K) . 

We shall need further the function u, which is defined modulo an additive con- 
stant by the relations: 

du du^ 
dxi dx2 ' 
du _ du0 
8x2 dxi 

The level sets of u form a congruence of curves in ft. The part of the boundary T j 
belongs to this congruence. We define the following system of open neighbourhoods 
named V{u), with the aid of this congruence: 

WA G V{u) dA \ Tu it is locally a level set of u . 

For any A eV{u) we denote by duA the part of the boundary of A belonging to 
or r^, i.e. 

duA = dAn {Tlurl) . 

The remaining part of dA is denoted by dfA. 

Let K be a rectifiable curve and CI' G V{u) such that 

K\rf c n'udu^' . 

For the couple {K, Q,') we introduce the following stress field: 




V'U0 in 

in n' . 
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This stress field is statically admissible with respect to the body with reference con- 
figuration 0,\K and boundary displacement no(l). Therefore we have the following 
inequality: 

J / I Vn/^ P dx > / (an) ■ uo{l) dH^ j \ a \^ dx . 

The latter inequality can be put in terms of Mumford-Shah functional / like this: 

i{(b)-i{K)< \ [ |Vn0|2 dx - GH\K\Tf) . 
^ Jn' 

The reason for which we have put 7i^{K\T f) instead of 7i^{K) is that uk = UK\rf 
but I{K) > I{K \ T f). As a consequence, sets K with a part on are disqualified 
to be minimizers of /. 

Let us denote by r the tangent vector field in direct sense to d^' . After few 
calculations we obtain from the previous inequality the estimation: 

1(0) - i{K) <\ f U0- (Vur) - GTiHK \Tf) . (6.2.4) 
^ Jan' 

We deduce that for any D > we have: 

£)2 I- 

i{^,D)-I{K,D) <— U0-{Vut)- Gn\K\rf) . (6.2.5) 
2 Jan' 

Let us return to the congruence of curves defined by u and consider the projection 
function on with respect to the congruence. For any set i? C F^ we denote by 
Vu{B) the variation of u on B. We see that: 

^ / (Vnen) • dH^ > ]-D^ Vu{P{K)) , (6.2.6) 
because of the equality: 

inf|— / (Vngn) •'U0 dW^ : rj' G y(u),K\Fj C Jl' j = --D^ yu(P(if)) . 
From the inequalities (|6.2.5jl and (|6.2.6jl we obtain the improved estimation: 

i{^,D)-I{K,D) < ^D^ Vu{P{K)) - GnHK\Tf) . (6.2.7) 
Let us remark that 

Vu{K) > Vu{P{K)) 

therefore we have: 

i{%,D)-i{K,D) < ^D'^VuiK) - GnHKXVf) . 
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Due to the assumption (recall the notation u{D) = u{uq{D),(!>)) 

C >| Vn(l)n p> c> (6.2.. 
we have Vu{K) < CH^{K) hence (if we suppose that K nVj is neghgible): 



i{9,D)-i{K,D) < 



-D^C -G 



(6.2.C 



Therefore, if 



-D'^C < G 
2 



from 1)6.2. 9|) we see that minimizes I{-,D), which proves the point i) of the theorem. 

Let us go back to ()6.2.7j) and introduce T{K) as the curve with the properties: 
pi) for any Q' G V{u), iiK\Tf cQ' then T{K) C n', 
p2) T{K) is a length minimizer in the class of curves that fulfills 

Pl). 

We remark that T{K) might not be unique, but it always exists. 
It is straightforward that T{K) = T{P{K)) and Vu{P{K)) = Vu{T{K)). We 
have then: 



i{%,D)-i{K,D) < ^^D"^ Vu{T{K)) - Gn\T{K)) 



From the assumption (|6.2.8|) we see that 

Vu{P{K)) - Gn^T{P{K)) > 



G 



nHviK)) 



(6.2.10) 



(6.2.11) 



Therefore, if 



-c > G 



then the right member of (|6.2.11() is positive and it attains the maximum when 
TC^{r{K)) is maximal. This happens when T{K) separates from F^. In this case 
is easy to see that we have equality in the relation 1)6.2. 1U() . which proves the point 
ii) of the theorem. 

The proof of iii) it is now straightforward. If C = c then 
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D'C 



7 Conclusions and perspectives 

The first model contains only a constant connected to fracture, namely the con- 
stant of Griffith G. The main qualities of this model are: 

i) crack appearance is allowed, together with crack propagation. 
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ii) there is no restriction concerning the pattern of the crack during 

its evolution. 

We have seen that in the first model the critical stress which lead to fracture (or 
crack appearance) is not a constant of material. 

The second model contains two constants of material connected to fracture: G 
and a constant with the dimension of a stress named S. In this model the critical 
stress which lead to crack appearance is a constant of material, related to S. This 
model has the same qualities as the first. 

These two models are fully macroscopical, in the sense that no fracture mecha- 
nism based on micro-cracks or other micro-defects was supposed. 

The main open theoretical problem is the general existence of an energy min- 
imizing movement according to our definitions. Below is described an existence 
result based on a sound physical assumption (|7.().lj) . Nevertheless, we do not know 
if HT.U.lf) can be proved from the basic assumptions of the model. 

Theorem 7.1. Let us consider for a given s an incremental solution k i— > 
{u^{k),K'^{k)) £ M, according to definition 4-1-- For any k € N we introduce the 
displacement 

v%k + l) = u{{k + l)/s,K'{k)) . 

Let us suppose that the power communicated by the rest of the universe to the body 
is bounded at any moment t. The incremental form of this assumption consists in 
the existence of a constant P such that for any k and s we have 

{TiK^k))^iuoiik + l)/s) + uoik/s)) ,uo{ik + l)/s)-uoik/s)) < P/s . (7.0.1) 

Then for any t > there exist diverging sequences (sj)j and {ki)i such that ki/si 
converges to t and {u^' , K^^){ki) converges to an element of M {u,K){t). 

Proof: Prom the minimality assumption on the incremental solution we have 
for any k £ N the inequality: 

J{iu'{k),K'{k)),{^,'{k + l),K'{k))) > J{{u'{k),K^{k)),{u'{k + l),K'{k + l))) . 
This inequality means that: 

[ w{Vv'{k + l)) dx > I wiVu'ik + l)) dx + Gn'^"\K'{k + l)\K'{k)) . 
Jn Jn 

The crack growth condition K^{k) C K^[k + 1) implies that the latter relation can 
be put in the following form: 

( / w{Vw'{k + l)) dx - [ w{Vu'{k)) dx] + (7.0.2) 
\Jn Jn J 

+ / w{Vu'{k)) dx + G7i'^'\K'{k) > I w(Vu'{k+l)) dx + GT-C'-HK'ik+l)) . 
Jn Jn 
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This is the incremental form of the Griffith criterion of crack propagation ()3.3.6|) . 
Indeed, we have the chain of equahties: 

[ w{Vv'{k + l)) dx - [ w{Vu'{k))dx = 
Jn Jn 

= ^{T{K%k))uo{{k + l)/s),uo{{k + l)/s)) - ^{T{K%k))uo{k/s),uo{k/s)) = 

= {TiK%k))^iuoiik + l)/s)+uoik/s)),uoiik + l)/s)-uoik/s)) . 

v'^(k + 1) represents the displacement of the body with the boundary displacement 
uo(A:/s + 1/s) in the presence of the crack K'^{k). u^(fc) represents the displacement 
of the body with the boundary displacement UQ^k/s) in the presence of the same 
crack K^{k). According to (|3.3.5() . the quantity 

j w{V^'{k + l)) - j wiVu^k)) dx^/Q 

is the discretized expression of the power communicated by the rest of the universe to 
the body at the moment k/s, when a time discretization with step 1/s is considered. 
We deduce from the inequality ()7.U.2() that 

P/s+ I w{Vu'{k)) dx + GW'HK'ik) > [ w{Vu%k+l)) dx + Gn''-\K'{k+l)) . 
Jn Jn 

We have therefore: 

Pk/s > [ w(yu'{k + i)) dx + Grr-^{K\k + i)) . 

Jn 

Prom the compactness theorem for SBD space and the latter inequality we deduce 
that for any t > Q there exist diverging sequences {si)i and {ki)i such that ki/si 
converges to t and K^^){ki) converges to an element of M {M^K){t). □ 

In the paper [AB] Ambrosio & Braides introduce a generalized minimizing move- 
ment based model for the propagation of a crack in the presence of viscous forces 
in the body. They give as initial datum at i = the anti-plane displacement 
uo G SBV(J7, i?) n L°^(Q,R). For a given s they recursively define a sequence 
in SBV(r2, ii) and an increasing sequence of closed rectifiable sets {Kf,)k as 
follows: Uq = no, Kq = and = -f^l+i = U K^, where w \s a minimizer 
of the functional 

j \ Vv\^ dx + rt-^{Sv\K'k) + s [ Iv-ull"^ dx (7.0.3) 
Jn Jn 

over the set of all v such that: 

V G SBV(i7,i?) , ||f||oo < ||no||oo • 
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The generalized minimizing movements obtained as limits of such incremental solu- 
tions, when s diverges, correspond to the following situation: a body evolves from 
the initial state uq, with the initial crack S^o, under a constant imposed boundary 
displacement. The equation of evolution for the displacement is: 

div Vu{t) + u{t) = . 

The authors obtain an existence result for the generalized minimizing movement 
introduced by them. After the introduction of the piecewise constant function: 

they find the following estimation: 

\\u'{t') -u'(t)\\L2 < M^t' -t+^ ift'>t . (7.0.4) 

Therefore there exists a diverging sequence {si)i such that u^"- converges to u uni- 
formly in L°^([0, T], L^(r2, i?)), for all T > 0. A consequence of this result is that 
the crack appearance is forbidden in this model. 

This result is obtained under the assumption of constant imposed boundary 
displacement, equal to the trace on the boundary of the initial datum uq. 

It is natural to introduce the Lame constant ^ and the viscosity A in the expres- 
sion of the functional ()7.0.3|) and modify it like this: 

[ filVvl"^ dx + W"~^(S^, \K|) + As / Iv-ull"^ dx . 
Jn Jn 

We obtain the more physical case of an anti-plane displacement satisfying at any 
moment t the equation: 

div iiVu{t) + \u{t) = . 
The estimation ()7.0.4() becomes 

\\u'{t')-uHt)\\L2 < M\ t'-t + ^ ift'>t . 

V As 

We expect to obtain our first model, in the case of anti-plane displacements, when 
the viscosity A converges to 0. It is easy to see that if A converges to then the 
uniform estimation from above is lost, hence there is no contradiction between the 
fact that in our model crack appearance is allowed and the fact that in the model 
of Ambrosio & Br aides crack appearance is forbidden. 

As a conclusion, an open direction of research consists in the use of more general 
minimizing movements in order to study the propagation of a crack in the presence 
of viscous effects (as is the paper [AB]) or in the case of an elasto-plastic body. 
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The models presented in the paper are of apphcative interest. In order to use 
them we have to know how to minimize a Mumford-Shah functional. This can 
be done by approximating, in the sense of variational convergence, the original 
functional by a less strange one. The idea is to replace the pair displacement- 
crack (u, K) with the pair (u, /), where / is a smoothed version of the characteristic 
function of the crack set K, taking values in the interval [0, 1] . The original functional 
may be replaced by an Ambrosio-Tortorclli approximation, introduced in [ATI], 
[AT2]. This opens the path to future interesting numerical results. 
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